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ABSTRACT 

Several independent lines of evidence now point to a connection between the physical 
processes that govern radio (i.e. jet) and X-ray emission from accreting X-ray binaries. 
We present a comprehensive study of (quasi-)simultaneous radioX-ray observations 
of stellar black hole binaries during the spectrally hard X-ray state, finding evidence 
for a strong correlation between these two bands over more than three orders of mag- 
nitude in X-ray luminosity. The correlation extends from the quiescent regime up to 
close to the soft state transition, where radio emission starts to decline, sometimes be- 
low detectable levels, probably corresponding to the physical disappearance of the jet. 
The X-ray transient V 404 Cygni is found to display the same functional relationship 
already reported for GX 339-4 between radio and X-ray flux, namely S ra dio oc S^ ' ■ 
In fact the data for all low/hard state black holes is consistent with a universal rela- 
tion between the radio and X-ray luminosity of the form L ra dio oc L x ' 7 . Under the 
hypothesis of common physics driving the disc-jet coupling in different sources, the 
observed spread to the best-fit relation can be interpreted in terms of a distribution 
in Doppler factors and hence used to constrain the bulk Lorentz factors of both the 
radio and X-ray emitting regions. Monte Carlo simulations show that, assuming little 
or no X-ray beaming, the measured scatter in radio power is consistent with Lorentz 
factors ^ 2 for the outflows in the low/hard state, significantly less relativistic than 
the jets associated with X-ray transients. When combined radio and X-ray beaming 
is considered, the range of possible jet bulk velocities significantly broadens, allowing 
highly relativistic outflows, but implying therefore severe X-ray selection effects. If the 
radio luminosity scales as the total jet power raised to x > 0.7, then there exists an 
X-ray luminosity below which most of the accretion power will be channelled into the 
jet, rather than in the X-rays. For x — 1.4, as in several optically thick jet models, the 
power output of 'quiescent' black holes may be jet-dominated below Lx - 4x 10~ 5 

LEdd ■ 
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1 INTRODUCTION 

There is strong observational evidence for the fact that pow- 
erful radio-emitting outflows form a key part of the accretion 
behaviour in some states of X-ray binary systems. Due to its 
high brightness temperature, 'nonthermal' spectrum and, in 
some cases, high degree of polarisation, radio emission from 
black hole binaries is believed to originate in synchrotron 

* egallo@science.uva.nl 



radiation from relativistic electrons ejected by the system 
with large bulk velocities (Hjellming & Han 1995; Mirabel 
& Rodriguez 1999; Fender, 2000, 2001a,b,c). 
Black hole binary systems are traditionally classified by their 
X-ray features (see Nowak 1995; Poutanen 1998; Done 2001; 
Merloni 2002 for recent reviews), namely: a) the relative 
strength of a soft 'black body' component around 1 keV, 
b) the spectral hardness at higher energies c) X-ray lumi- 
nosity and d) timing properties. Different radio behaviour 
is associated with several 'X-ray states', according to the 



2 E. Gallo, R. P. Fender & G. G. Pooley 



following broad scheme. The low/hard state is dominated 
by a power-law spectrum, with a relatively low luminos- 
ity and an exponential cut-off above about 100 keV and 
little or no evidence for a soft, thermal component. It is 
associated with a steady, self-absorbed outflow that emits 
synchrotron radiation in the radio (and probably infrared) 
spectrum. The quiescent/off state, characterised by an ex- 
tremely low X-ray flux, may simply be interpreted as the 
hard state 'turned down' to lower accretion rates and radia- 
tive efficiency. X-ray spectra from high/soft state Black Hole 
Candidates (BHCs) are instead dominated by thermal radi- 
ation, while the core radio emission drops below detectable 
levels, probably corresponding to the physical suppression 
of the jet. In the very high state both the thermal and the 
power law components contribute substantially to the spec- 
tral energy distribution. At a lower luminosity level an in- 
termediate state is also observed, with properties similar to 
those of the very high state. For both the very high and the 
intermediate state the connection with radio behaviour is 
not yet clearly established. Corbel et al. (2001) show that 
the radio emission from XTE 1550-564 in the intermediate 
state was suppressed by a factor > 50 with respect to the 
hard state, while Homan et al. (2001) claim that intermedi- 
ate and very high states can actually occur at a wide range 
of luminosities. 

Transitions between states are often associated with multi- 
ple ejections of synchrotron emitting material, possibly with 
high bulk Lorentz factors (Hjellming & Han 1995; Kuulkers 
et al. 1999; Fender & Kuulkers 2001). 

As already mentioned, BHCs in the low/hard state, like 
Cygnus X-l and GX 339-4, are characterised by a flat or 
slightly inverted radio spectrum (a = A log S v /A log v ~ 0), 
interpreted as arising from a collimated, self-absorbed com- 
pact jet, in analogy to those observed in active galactic nu- 
clei (Blandford & Konigl 1979). With the direct imaging of 
a resolved compact radio jet from the core of Cygnus X-l 
(Stirling et al. 2001), this association has been confirmed. 
Radio emission from X-ray binaries, especially the BHCs, is 
increasingly interpreted as the radiative signature of jet-like 
outflows. 

It has been generally accepted that the soft thermal com- 
ponent of BH spectra originates in an optically thick, ge- 
ometrically thin accretion disc (Shakura & Sunyaev 1973), 
whereas the power law component is produced by Comptoni- 
sation of 'seed' photons in a hot, rarefied 'corona' of (quasi-) 
thermal electrons (Shapiro, Lightman & Eardley 1976; Sun- 
yaev & Titarchuk 1980; Haardt & Maraschi 1991; Poutanen 
& Svensson 1996). Although this picture can successfully 
reproduce the X-ray behaviour, it can not yet address the 
clear correlation between radio and X-ray emission estab- 
lished for the persistent BHCs GX 339-4 and Cygnus X-l 
while in the hard state (Hannikainen et al. 1998; Brocksopp 
et al. 1999; Corbel et al. 2000; Corbel et al. 2003). More- 
over, some hard state sources, like XTE Jl 118+480 and GX 
339-4, show evidence for a turnover in the infrared-optical 
band, where the flat -to-inverted radio spectrum seems to 
connect to an optically thin component extending up to the 
X-rays (Corbel & Fender 2002; Markoff et al. 2003a,b and 
references therein), suggesting again that the jet plays a role 
at higher frequencies. 

Hence, all the evidence points to the corona in these sys- 
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Figure 1. Radio against X-ray flux density, scaled to a distance 
of 1 kpc and absorption corrected, for V 404 Cygni and GX 339- 
4. Lines denote the fits to the datasets: short and long dashed for 
GX 339-4 and V 404 Cygni respectively. It is found that the data 
of V 404 Cygni are well fitted by the same functional relationship 
reported by Corbel et al. (2003) for the BHC GX 339-4, that is 

S ra dio OC (S X ) - 7 . 



tems being physically related to the presence of a jet: by far 
the simplest interpretation therefore is that the Comptonis- 
ing region is just the base of the relativistic outflow (Fender 
et al. 1999b; Merloni & Fabian 2002; Markoff et al. 2003a). 
However, joining these two previously independent scenarios 
is somewhat problematic because they often require differ- 
ent electron distributions and geometries. 
Due to the fast timescales in X-ray binary systems, only 
simultaneous radio and X-ray observations provide the nec- 
essary tools to probe this conjecture. The following results 
extend and complete those presented in Gallo, Fender & 
Pooley (2002). 



2 THE SAMPLE 

Our aim was to compile (quasi-)simultaneous radio and X- 
ray observations of BHCs during the low/hard state. To 
this purpose, we have collected all the available (to our 
knowlege) data from the literature and made use of our 
own simultaneous observations as well. These were taken 
with the Ryle telescope at 15 GHz (see Pooley & Fender 
1997 for more details) and combined with one-day aver- 
ages from RXTE ASM (this refers to Cygnus X-l, Cygnus 
X-3 and GRS 1915+105). Table 1 lists the information for 
the ten low/hard state BHCs for which we have at our dis- 
posal (quasi-)simultaneous radio and X-ray coverage (see 
Section 5, Table 3 for 'non-canonical' hard state sources, 
such as Cygnus X-3 and GRS 1915+105): distance, mass, 
orbital inclination, measured hydrogen column densities (see 
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Table 1. System parameters for the ten hard state BHCs under consideration. Distance and Njj references are given in parentheses next 
to each value. Inclination and BH mass estimates, unless differently specified, are all taken from Orosz (2002). The last column refers to 
the literatures' sources from which we have obtained (quasi-) simultaneous radio and X-ray fluxes; no reference appears in case of our 
own observations (Cygnus X-l). 



Source 


Distance(ref) 


Inclination(12) 


BH mass(12) 


Nff(ref) 


Data 




(kpc) 


(degree) 


(Mo) 


(10 21 cm- 2 ) 




Cygnus X-l 


2.1 (1) 


35±5 


6.85-13.25 


6.2 (16) 




V 404 Cygni 


3.5 (2) 


56±4 


10.06-13.38 


5.0 (17) 


21,22 


GRS 1758-258 


8.5 (3) 


? 


~8-9 (14) 


14.0 (3) 


23 


XTE J1118+480 


1.8 (4) 


81±2 


6.48-7.19 


0.1 (18) 


24,25 


GRO J0422+32 


2.4 (5) 


44±2 


3.66-4.97 


2.0 (5) 


24 


GX 339-4 


4.0 (6) 


15-60 (13) 


5.8±0.5 (15) 


6.0 (6) 


26 


IE 1740.7-2942 


8.5 (7) 


? 


? 


118 (19) 


27 


XTE J1550-564 


4.0 (8,9*) 


72±5 


8.36-10.76 


8.5 (20) 


20,28 


GS 1354-64 


10.0 (10) 


? 


? 


32.0 (9) 


29 


4U 1543-47 


9.0 (11) 


20.7±1.5 


8.45-10.39 


3.5 (11) 


30 



References : 1: Massey et al. 1995; 2: Zycki, Done & Smith 1999; 3: Main et al. 1999; 4: McClintock et al. 2001; 5: Shrader et 
al. 1997; 6: Zdziarski et al. 1998; 7: Sunyaev et al. 1991; 8: Kong et al. 2002; 9: Tomsick et al. 2001; 10: Kitamoto et al. 1990; 11: 
Orosz et al. 1998; 12: Orosz 2002; 13: Cowley et al. 2002; 14: Keck et al. 2001; 15: Hynes et al. 2003; 16: Schulz et al. 2002; 17: 
Wagner et al. 1994; 18: Dubus et al. 2001; 19: Gallo & Fender 2002; 20: Tomsick et al. 2001; 21: Han & Hjellming 1992; 22: 
Hjellming et al. 2000; 23: Lin et al. 2000; 24: Brocksopp et al. 2003; 25: Markoff, Falkcc & Fender 2001; 26: Corbel et al. 2000; 27: 
Heindl, Prince & Grunsfcld 1994; 28: Corbel et al. 2001; 29: Brocksopp et al. 2001; 30: Brocksopp, private communication. 
*For XTE J1550-564 a distance of 4 kpc is assumed by both Kong et al. (2002) and Tomsick et al. (2001), as average value between 
2.5 and 6 kpc, given by Sanchez-Fernandez et al. (1999) and Sobczak et al. (1999) respectively. 



Section 2.1) and literature references. 

Both the X-ray and the radio intensity values come from 
several different instruments and telescopes. Radio flux den- 
sities have been measured in different frequency bands, rang- 
ing from 4.9 up to 15 GHz; nevertheless we generically refer 
to 'radio flux densities' based on the evidence that, while 
in the low/hard state, black hole radio spectra are charac- 
terised by almost flat spectra (a ~ 0) spectral index (Fender 
2001a). 

X-ray fluxes, taken either from spectral fits or from light 
curves, have been converted into Crab units in order to 
be easly compared with radio flux density units (1 Crab ~ 
1060 /iJy; energy range 2-11 keV). For this purpose, X-ray 
fluxes/luminosities in a given range have been first converted 
into corresponding values between 2-11 keV, and then ex- 
pressed as flux density. For those sources whose X-ray flux 
has been derived from count rates, the conversion into Crab 
has been performed according to the factors provided by 
Brocksopp, Bandyopadhyay & Fender (2003). 



2.1 Absorption corrections 

Whenever X-ray flux density has been evaluated from count 
rates or absorbed fluxes, we wanted to compensate for ab- 
sorption by calculating the ratio between the predicted flux 
from a hard state BH with a measured Nh value, and the 
predicted flux corresponding to no absorption, as follows. 
We have first simulated with XSPEC typical spectra of hard 
state BHCs as observed by Chandra ACIS for ten diffent val- 
ues of hydrogen column density ranging from zero up to 12.5 
xlO cm~ . A 'typical' spectrally hard BH's spectrum is 
well fitted by an absorbed power law with photon index 1.5. 
By keeping fixed the flux corresponding to no absorption, 
the points turn out to be well fitted by a simple exponential 



relation, which allows to express the ratio F a bs/F una b s as 
follows: 



Fab 



Funab 



= exp 



-(Nh/W 22 cm- 2 ) 



18.38 



(1) 



The procedure described has been applied to X-ray fluxes 
below the transition luminosity between hard and soft state 
(see Section 3.2). Above that value, the spectrum is not re- 
produced by a simple power law. In this regime, the X-ray 
spectrum is usually well fitted by an absorbed power law 
with photon index T ~ 2.4 plus a disc blackbody compo- 
nent, with a typical temperature of around 1 keV. Since in 
this case the 2-11 keV spectrum is almost entirely domi- 
nated by thermal emission, the previous simulations have 
been repeated for soft state BHCs by approximating the 
spectrum with a disc blackbody emission at 1 keV. We have 
obtained: 



= exp 



-{N H /10 22 cm- 2 ) 



8.67 



(2) 



The latter correction has been applied to detections above 
the hard-to-soft state transition. 



3 RADIO VS. X-RAY FLUX DENSITIES 
3.1 GX 339-4 and V 404 Cygni 

In Figure 1 we plot radio against X-ray flux densities (mjy 
vs. Crab), scaled to a distance of 1 kpc and absorption cor- 
rected, for GX 339-4 and V 404 Cygni, the two sources for 
which we have at our disposal the widest coverage in terms 
of X-ray luminosity. 

GX 339-4 was discovered as a radio source by Sood & 
Campbell-Wilson (1994). When in the low/hard state, it 



4 E. Gallo, R. P. Fender & G. G. Pooley 

10J IP" 5 

I I I I l l l | I I I I I l l l | I l 

A V404 Cyg (3.5 kpc) 

10 3 h- • GX 339 4 ( 4 k P c ) 

# XTE ,11 1 18 + 480 (1.8 kpc) 
A 4U 1543-47 
x GS 1354 64 (10 kpc) 
100 h # IE 1740.7 2942 (8.5 kpc) 

Ci XTE J 1 550 564 (4 kpc) 



■ 

W 

CD 

3 



+ 

CJ 



10 



10 4 0.001 

I I I I I I I I I I I I 

L/L Kdd (10 MJ 



0.01 



m 



m 




u CRO J0422+32 

_ GRS 1758-258 




■ 




S .. =kx(S Y ) +0; 

radio v A- 1 

k=223±156 



0.00 I 



I I I I I 

0.1 



1 
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Figure 2. Radio flux density (mjy) is plotted against X— ray flux density (Crab) for a sample of ten hard state BHs (see Table 1), scaled 
to a distance of 1 kpc and absorption corrected (this means that the axes are proportional to luminosities). On the top horizontal axis 
we indicate luminosity, in Eddington units for a IOMq BH, corresponding to the underlying X— ray flux density. An evident correlation 
between these two bands appears and holds over more than three orders of magnitude in luminosity. The dashed line indicates the best-fit 
to the correlation, that is S ra di = k x (Sx) + ' , with k = 223 ± 156 (obtained by fixing the slope at +0.7, as found individually for 
both GX 339-4 and V 404 Cygni; see Section 4.1). Errors are given at 3— a confidence level, and arrows also represent 3— a upper limits. 



is characterised by a flat or slightly inverted (a it 0) ra- 
dio spectrum (see Corbel et al. 2000) and its synchrotron 
power has been shown to correlate with soft and hard X- 
ray fluxes (Hannikainen et al. 1998; Corbel et al. 2000). By 
means of simultaneous radio:X-ray obsevations of GX 339- 
4, Corbel et al. (2003) have recently found extremely inter- 
esting correlations between these two bands: in particular, 
Ss.sghz oc S 3 t_Qkly' 01 (slightly different slopes - within the 
hard state - have been found depending on the X-ray en- 
ergy interval). When fitted in mjy vs. Crab (scaled to 1 kpc 
and absorption corrected), the relation displays the form: 

Sradio = kaX339-4 X (Sx) + 0i71±0 - 01 (3) 

kGX339-4 = 126 ± 3 

The correlation appears to hold over a period of three years 
- 1997 and between 1999-2000 - during which the source 
remained almost constantly in a spectral hard state (with a 



transition to the high/soft state, Belloni et al. 1999, when 
the radio emission declined below detectable levels). Fig- 
ure 1 shows radio against X-ray flux densities of GX 339-4 
corresponding to simultaneous ATCA/RXTE observations 
performed between 1997 and 2000 (Corbel et at 2000, 2003). 
Note that points above 1 Crab (scaled), which all correspond 
to RXTE-ASM detections, clearly show a sharp decreas- 
ing in the radio power (see next Section). The correlation 
reported by Corbel et al. (2003) actually refers to RXTE- 
PCA data only; it is worth mentioning that, when ASM 
detections below 1 Crab (i.e. below the radio quenching) 
are fitted together with PCA points, the final result is con- 
sistent, within the errors, with the fit reported by Corbel on 
PCA data alone (that is, a slope of 0.70+0.06 is obtained in 
this case). 

Remarkably, we have found that detections of V 404 Cygni, 
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the source for which we have at our disposal the widest 
radio:X-ray coverage, are well fitted by the same functional 
relationship - albeit with no apparent cutoff - as GX 339-4 
(see Figure f). 

V 404 Cygni belongs to the class of X-ray transients, sources 
undergoing brief episodic outbursts during which their lumi- 
nosity can increase by a factor ~ 10 6 compared to periods 
of relative quiescence. All V 404 Cygni data - except for 
the quiescent lowest point - come from simultaneous radio 
(VLA; Han & Hjellming 1992) and X-ray {Gmga; Kitamoto 
et al. 1990) observations during the decay following its May 
1989 outburst, during which the source, despite very high 
and apparently saturated luminosity, never entered a spec- 
tral soft state and always maintained a very hard X-ray 
spectrum (Zycki, Done & Smith 1999). According to Hjellm- 
ing et al. (2000), the quiescent state of V 404 Cygni, since 
it ended the long decay after its 1989 outburst, has been 
associated with a 0.4 mjy radio source 1 . Quiescent X-ray 
flux refers to a 1992 measurement (Wagner et al. 1994 re- 
port 0.024±0.001 count/sec with ROSAT-PSPC), i.e. well 
before the onset of significant X-ray variability (see Kong et 
al. 2002 for details). 

Denoting S ra dio as the radio flux density in mjy and Sx as 
the X-ray flux density in Crab, we have obtained: 

Sradio = kv404 X (Sx) + 0,70±() ' 20 (4) 

AV404 = 301 ± 43 

The Spearman's rank correlation coefficient is 0.91; the 
two sided significance of its deviation from zero equals 4.2 
xl0~ 3 . 

These results indicate that S ra dio oc (Sx)"" ' 7 is a funda- 
mental property of the radio:X-ray coupling in the hard 
state, rather than a peculiarity of GX 339-4. ft is worth 
stressing that the fitted slopes for V 404 Cygni and GX 
339-4 are identical within the errors, with the same nor- 
malisations within a factor 2.5, while detections from other 
sources below 1 Crab (scaled), although much narrower lu- 
minosity ranges, are all consistent with the same placing in 
the radio:X-ray plane, as discussed in the next Section. 



3.2 Broad properties of the correlation 

In Figure 2 we plot radio flux densities (mjy) against X-ray 
fluxes (Crab), scaled to a distance of 1 kpc and absorption 
corrected, for all the ten hard state BHCs listed in Table 1. 
Note that this scaling and correction means that the axes 
are proportional to luminosities. 

Besides GX 339-4, Cygnus X-l displays a positive S ra dio '■ 
Sx correlation followed by a radio turnover around 3 Crab, 
whereas three other sources, namely XTE J1118+480, 4U 
1543-47 and GS 1354-64, lie very close to the relations in- 
ferred for GX 339-4 and V 404 Cygni. The remaining four 
systems (GRS 1758-258, GRO J0422+32, IE 1740.7-2942 



1 Starting in early 1999, VLA obervations showed fluctuations 
ranging from 0.1 to 0.8 mjy on time scales of days; even more 
extreme radio fluctuations in February 2000 were accompaincd 
by strong variability in the X-ray band as well (Hjellming et 
al. 2000). 




Scaled X-ray flux density (Crab) 

Figure 3. Top panel: Ryle telescope at 15 GHz and RXTE ASM 
daily averages for Cygnus X-l between January 1996 and January 
2003. Fluxes have been scaled to a distance of 1 kpc. Bottom 
panel: open squares for HR (Hardness ratio)l, defined as the ratio 
between the count rates at 3-5 keV and 1.5-3 keV. Filled circles 
for HR2, defined as the ratio between 5-12 keV and 3-5 keV. 
A softening of the X-ray spectrum is shown to correspond to a 
quenching in the radio emission. 

and XTE J 1550-564) seem instead to have already under- 
gone the radio quenching. 

We can assert that these ten BH candidates display very 
similar behaviour in the S ra dio vs. Sx plane. There is evi- 
dence for a positive radio:X-ray correlation over more than 
three orders of magnitude in terms of Eddington luminosity, 
as indicated in the top horizontal axis, where we show the 
LjljEdd ratio corresponding to the underlying X-ray flux 
(Lsdd — 1.3 x 10 39 erg/sec for a 10 M BH). As an example, 
the total luminosity of Cygnus X-l in the 0.1 to 200 keV 
band, while the source is in the low/hard state, is ~ 2% 
of the Eddington luminosity for a 10 M a BH (di Salvo et 
al. 2001); that corresponds to mean flux of about 2.5 Crab 
(scaled to 1 kpc). 

3.3 Jet suppression in the soft state 

As already mentioned, approaching the soft state radio emis- 
sion from BHCs seems to be quenched below detectable lev- 
els. Such a behaviour is well recognizable in at least two 
sources of our sample. As visible in Figure 2, both Cygnus 
X-l and GX 339-4 display a clear turnover after which the 
radio power dramatically declines and reaches undetectable 
levels within a factor of two in X-ray flux (see also Fender et 
al. 1999b). In addition, XTE J1550-564, although suffering 
from poor statistics and large uncertainty in distance esti- 
mates (2.5-6 kpc), has been detected during the hard state 
(Corbel et al. 2001; Tomsick et al. 2001) at 'half-quenched' 
level (about 1 Crab:20 mjy scaled) while the radio emission 
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dropped down significantly in the intermediate/very high 
state (about 6 Crab:<0.15 mjy scaled). 
The values of the X-ray flux density (scaled) correspond- 
ing to the turnover in radio flux vary between the sources: 
GX 339-4, about 1 Crab; Cygnus X-f, about 3 Crab; V 
404 Cygni > 1 Crab. This difference could be related to 
differences between the parameters that govern the pow- 
ering/quenching mechanism(s) of the radio emitting jet. A 
discriminant parameter might be the BH mass, which is es- 
timated to be ~10 Mq in the case of Cygnus X-f while it 
has been recently constrained around 6 Mq in the case of 
GX 339-4 (Hynes et al. 2003). If this hypothesis is correct, 
we would expect not to see jet quenching at X-ray luminosi- 
ties below the Cygnus X-f turnover for those systems whose 
BH mass has been estimated to be bigger than 10 Mq, as, 
for instance, the case of V 404 Cygni, which does not ap- 
pear to be quenched in radio up to about 3 Crab, and is 
known to possess a f0-f4 M BH (Orosz 2002, Hjellming et 
al. 2000). Despite the poor statistics, this picture seems to 
bear out the hypothesis that the X-ray luminosity at the ra- 
dio quenching might positively correlate with the BH mass, 
being consistent with the jet suppression occurring at a con- 
stant fraction - a few percent - of the Eddington rate. This 
is explored further in Gallo, Migliari & Fender (in prepara- 
tion). 

In Figure 3 the quenching of the radio power in Cygnus X- 
f (monitored simultaneously in radio and X-ray between 
January f 996 and January 2003) as a function of the X-ray 
flux density is shown together with X-ray hardness ratios. 
The 3-f2 keV X-ray spectrum softens until about 5 Crab 
(scaled), whereas radio quenching begins around 3 Crab 
(scaled). While it is clear that the quenching occurs some- 
where near the point of transition from low/hard to softer 
(intermediate or high/soft) states, pointed observations will 
be required to see exactly what is happening to the X-ray 
spectrum at the quenching point. 

4 SPREAD TO THE CORRELATION 
4.1 Best— fitting 

So far we have established two main points by looking at 
the distribution of low/hard state BH binaries in the radio 
vs. X-ray flux density plane: 

• Independently of the physical interpretation, S ra dio fx 
(Sx) +0 - 7 for GX 339-4 and V 404 Cygni from quiescence up 
to close to the hard-to-soft state transition. All other hard 
state BHCs lie very close to these correlations, with similar 
normalisations. 

• At a luminosity of a few percent of the Eddington rate, 
close to the hard-to-soft state transition, a sharp turnover 
is observed in the radio:X-ray relation, that is, the radio 
flux density drops below detectable levels. 

Bearing this in mind, our purpose is now trying to find a re- 
liable expression for a 'best-fit' relationship to all hard state 
BHs and to estimate the spread relative to such a relation. 
Assuming 0.7 as a universal slope during the low/hard state, 
we have determined the normalisation factor by fitting all 
the data - Cygnus X-f excluded (see comments at the end 



of this Section) - up to f Crab, below which quenching does 
not occur for any system. 

In this way we are able to provide an empirical relationship 
valid for all the hard state BHs, that we will call 'best-fit' 
in the following. We have obtained: 

Sra^o = k x (Sx) +0 - 7 , fc = 223 ±156 (5) 

The best-fit and its spread are indicated in Figure 2 in 
dashed and dot-dashed lines respectively. A scatter of about 
one order of magnitude in radio power is particularly inter- 
esting, especially in view of comparing the observed spread 
to the one we expect based on beaming effects (see next Sec- 
tion). 

The choice of excluding the whole dataset of Cygnus X-f 
is related to its unusual behaviour in the radio:X-ray flux 
plane. In fact, a visual inspection of Figure 2 already sug- 
gested that the points below the radio power quenching be- 
long to a line with a steeper slope than 0.7. In addition, 
despite its relatively low inclination to the line of sight (see 
next Section for clarity), Cygnus X-f lies on the lower side 
of the correlation. It is possible that these characteristics 
can be explained in terms of strong wind absorption. It has 
been demonstrated that the wind from the donor OB star 
in Cygnus X-f partially absorbs the radio emission, up to 
about 10% (Pooley, Fender & Brocksopp, f999; Brocksopp, 
Fender & Pooley, 2002). In addition, since the jet bulk ve- 
locity during hard state is likely to be relatively low (see 
Section 4.3), approaching and receding jets contribute a 
similar amount to the total radio luminosity. However, in 
Cygnus X-l only a one-sided jet has been detected (Stirling 
et al. 2001); this means that, possibly, a significant fraction 
of the receding jet is lost through wind absorption (however, 
since the jet structure is about f 00 times bigger than the or- 
bit, it is unlikely that the wind could absorb the entire power 
of the receding component). Furthermore, because the flat- 
spectrum radio emission (corresponding to steady jet and 
generally associated with hard X-ray spectrum) is optically 
thick, this implies that when the jet power decreases, its size 
might also decreases linearly with flux. As the jet in Cygnus 
X— 1 becomes smaller, the S ra dio'-Sx relation will be subject 
to increasingly strong wind absorption, with the net result 
of steepening the correlation. 



4.2 Constraining the Doppler factors 

What kind of physical information can we deduce from the 
observed relation? 

Let us assume a very simple model, in which the same 
physics and jet/corona coupling hold for all hard state BHs 
and the observed functional relationship is intrinsic; then, 
one would predict the following placing in the S ra dio vs. Sx 
parameter space: 

(i) All sources lying on a line with the same slope if nei- 
ther X-ray or radio emission were significantly beamed (i.e. 
low Doppler factors). 

(ii) Different sources lying on lines with the same slope 
but different normalisations if X-rays were isotropic while 
radio was beamed, with radio-brighter sources correspond- 
ing to higher Doppler factors. 
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(iii) As point (ii) but with higher Doppler factors sources 
being brighter in both radio and X-ray if both were beamed. 

Despite the relatively small sample, we are able to place 
some constraints on these possibilities. In the first two sce- 
narios, where X-rays are isotropically emitted while radio 
power is beamed, we can express the observed radio lu- 
minosity as the product of the intrinsic (rest-frame) radio 
power times the effective Doppler factor A ra£ jio, defined as 
a function of approaching and receding Doppler factors 2 . If 
S, 

write 



radio,intr = k x (Sx) +0 ' 7 for all hard state BHCs, we can 



+0.7 



(6) 



Sradio.obs ^radio X S ra dio,intr ^radio X N X (>Sx) 

Assuming the same coupling for all sources - that means 
same normalisation N - the ratio S ra dio,i/ S ra dio,2 between 
the observed radio powers from source 1 and 2, at a fixed 
X-ray luminosity, will correspond to the ratio between their 
relative effective Doppler factors. 

Returning to the case of V 404 Cygni and GX 339-4, 
where S ra dio,V404/ ' S ra dio,GX339-4 ~ 2.5, we are drawn to 
the conclusion that GX 339-4, whose inclination is poorly 
constrained between 15 < i < 60° (Cowley et al. 2002), 
is likely to be located at a higher inclination than V 
404 Cyg, well established to lie at 56 ± 4° (Shahbaz et 
al. 1994). In fact, assuming that Fgx339-4 — IV404, a ra- 
tio Sradio,v404/Sradio,GX339~4 > 1 can be achieved only if 

«GX339-4 > 2V404- 

Clearly, the previous arguments are based on the assump- 
tion that the binary system inclination to the line of sight 
coincides with the inclination of the jet, while recent find- 
ings (Maccarone 2002) show that this is not always the case 
(the misalignment of the disc and the jet has been already 
observed in Galactic jet sources GRO J1655-40 and SAX 
J1819-2525). 

A recent work by Hynes et al. (2003) shows dynamical ev- 
idence for GX 339-4 being a binary system hosting a BH 
with mass 5.8 ± 0.5 M©. Interestingly, based on the spec- 
troscopic analysis by Cowley et al. (2002), see their Figure 
8, this value is also consistent with a high inclination of the 
system to the line of sight. 



4-2.1 Monte Carlo simulations I: radio beaming 

In order to link the measured scatter in the radio:X-ray re- 
lation to the beaming effects, and possibly constraining the 
Doppler factors, we have performed Monte Carlo simulations 
according to the following scheme. We have considered the 
four sources (see Table 2) whose radio:X-ray flux densities 
have been utilised to obtain the best-fit relationship below 
1 Crab; for each of them we have at our disposal a num- 



points n = Yl* =1 nj = 23. For each source j, a random incli- 
nation within the measured range has been generated and 
associated to an array whose dimension equals the number 

2 A radio := [(Sapp) 2 + (<5 rec ) 2 ]/2, where: 

Srec/app = T" 1 X (1 ± (3 COS 6) ~ 1 ; 

(i = v/c, is the bulk velocity of the radio— emitting material; 

r = (1 — /3 2 ) -0 - 6 ; 9 is the inclination respect to the line of sight. 



~fcf 



• random inclinations 
o measured inclinations 



A: = \5 Z +5 Z 1/2 

L app rec-l/ 



a(k)/k = 0.70 



CO 

d 



0.25 



0.5 



0.75 



Figure 4. Assuming a simple model in which radio emission is 
beamed while X-rays are isotropic, equation 6 implies that the 
spread in radio power to the best fit can be due to the distribution 
in Doppler factors. The simulated spread in A ra£ j io (defined as 
the ratio of standard deviation over the mean value), is plotted 
for 15 values of the jet bulk velocity between and 0.998c, and 
compared to the measured scatter in the radio:X-ray correlation. 
In order to maintain the spread in Doppler factors equal/smaller 
than observed, the outflow velocity must be smaller than 0.8 times 
the speed of the light, that is T ra dio < 2. Open circles represent 
Monte Carlo simulation run by allowing the inclination angles 
to vary within the observed ranges; filled circles are for random 
inclinations between 0-90° . 



Table 2. Observed sources for which mean Doppler factors have 
been simulated for 10 different values of the jet bulk velocity; mea- 
sured range of inclinations and number of simultaneous radioX- 
ray detections are listed. 



Source 


Inclination 


Number of 




(degree) 


detections (nj) 


V 404 Cygni 


56±4 a 


7 


XTE J1118+480 


81±2 a 


1* 


4U 1543-47 


20.7±1.5 a 


1 


GX 339-4 


15-60 6 


14 



* This single point corresponds to the average of 33 nearly 
simultaneous detections over a very narrow interval in X-ray 
luminosity. 

References: a: Orosz 2002; b: Cowley et al. 2002. 

of detections rij, for a total of 23 random inclinations. Then, 
for 15 different values of the outflow bulk velocity f3, the cor- 
responding Doppler factors Ap(j) have been calculated by 
running the simulation 10 4 times (i.e. for total of 23 x 10 4 
values for each f3) . After estimating the mean Doppler factor 
and its relative standard deviations (a (A)) for each value of 
[3, we have compared the simulated spread (defined as the 
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ratio <r(A)/A mea „ to the measured spread in normalisation, 
that is a(k)/k = 156/223 = 0.70. 

The result is shown in Figure 4: in order to keep the spread 
in Doppler factor equal/smaller 3 than observed (actually 
< 0.77 = 0.70+10%), the bulk velocity of the radio-emitting 
material must be lower than 0.8c, that is the Lorentz factor 
must be smaller than 1.7. 



For comparison, we ran the simulation allowing the incli- 
nation angles to vary between 0-90° (filled circles in Fig- 
ure 4); this is actually crucial in the light of what has 
been discussed by Maccarone (2002) about jet-disc misalign- 
ment in BH binaries, and also takes into account possible 
model-dependent errors in the estimation of i. The simu- 
lated spread in Doppler factors starts to significantly deviate 
from the previous one (calculated allowing the inclination 
angles for each source to vary within the measured values; 
open circles in Figure 4) only for very high bulk velocities. As 
expected, even if inclinations as and 90° (i.e. those incli- 
nations which translate into the highest and lowest possible 
values of cosine, respectively) are also taken into account, 
this strongly influences the mean Doppler factor only for 
really high values of /3. 

4-2.2 Monte Carlo simulations II: adding beamed X-rays 

Following a similar approach, it is possible to constrain the 
Doppler factor due to the combination of beamed radio and 
X-ray radiation. As before, we assume 0.7 as an intrinsic 
slope which relates radio and X-ray emission. Supposing 
that both X-rays and radio are beamed, we will write: 

Sradio, obs ^radio /~\ 

(Sx, ob sr 7 K (Ax) - 7 1 ' 

Note that the effective Doppler factor for the X-ray radia- 
tion will be defined as Ax = [(8 r ec) 2 ' 5 + (S apP ) 2 ' 5 ]/2 (assum- 
ing continuous ejection and photon index 1.5 for low/hard 
state; see Mirabel & Rodriguez 1994, equations 8, 9). 
Therefore, Monte Carlo simulations have been run by vary- 
ing both radio and X-ray bulk velocities with logarithmic 
steps in Lorentz factors, Tx and Y radio , between 1 and 30. 
The mean values of A ra di /(Ax) ' 7 together with their stan- 
dard deviations and spreads have been calculated for any 
combination of the two factors. The results are shown in 
Figure 5. 

Filled bold circles indicate those combinations of Y ra dio and 
Yx for which the simulated spread with inclinations varying 

3 Smaller than measured spread are also 'allowed' on the ground 
that errors in the distance estimates are likely to influence the 
observed distribution, causing an additional source of scatter to 
the relation, which is of course not related to any boosting effect. 
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Figure 5. Same scheme as Figure 4, but assuming a model in 
which both radio and X-ray beaming are allowed; in this case, 
Monte Carlo simulations have been performed with logarithmic 
steps in Lorentz factors (rather than linear in 0). The combina- 
tion of radio plus X-ray beaming has the effect of broadening the 
range of allowed Y ra dio with respect to the case of isotropic X- 
ray emission. In this picture open circles indicate combinations 
of Y ra dio and Yx for which the simulated spread with inclina- 
tions varying within the measured ranges is smaller than/equal 
to 0.77. Filled circles, instead, correspond to 'allowed' combina- 
tions of Y ra dio and Tx when the inclination to the line of sight is 
randomly chosen between 0—90° . In this case the range of possible 
combination is smaller due to the fact that extreme inclinations, 
such as or 90° arc also taken into account. 



between 0-90° is smaller than/equal to 0.77. Filled PLUS 
open circles, instead, correspond to 'allowed' combinations 
when the inclination to the line of sight varies within the 
measured ranges. 

In order to see why the beaming of both wavebands al- 
lows a wider spread of Y, consider the effect on the po- 
sition of a source in the flux-flux diagram: a Doppler- 
boosted, or deboosted, radio source moves parallel to the 
Sradio axis, while a source boosted or deboosted in both 
wavebands moves roughly (not exactly) parallel to the line 
Sradio = const x Sx , and therefore roughly along the di- 
rection of the observed correlation. Consequently beaming 
in both wavebands disturbs the correlation less. However, if 
X-rays were really highly beamed, that would imply strong 
X-ray selection effects in detecting BH binaries. 
Clearly, a possible independent estimation of the bulk veloc- 
ity of the X-ray emitting material (e.g. Beloborodov 1999; 
Maccarone 2003) would naturally allow a much narrower 
constraint on the Doppler factor of the radio-emitting ma- 
terial (see discussion for further details). 



These remarks are of course valid under the basic assump- 
tions that no beaming is affecting the X-rays. In addition 
we are considering a simple model in which both the bulk 
velocities and opening angles of the jets are constant (only 
under the latter assumption the probability of observing a 
source with a given inclination 9 is uniformly distributed in 
cos 9). 
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Figure 6. In addition to hard state BHs (open stars) we plot both single simultaneous radio:X-ray fluxes from black hole transients 
and detections from two canonical soft state plus somewhat 'extreme' sources, close to the Eddington accretion regime. Filled squares 
correspond to single outbursts from different sources (from Fender & Kuulkcrs 2001); filled circles are for persistent soft/high state 
sources LMC X— 1 and LMC X-3, upper limits; open polygons refer to GRS 1915+105, while open circles denote Cygnus X— 3 points. 
Based on the 0.7 correlation, and assuming that the radio luminosity scales as the total jet power raised to x = 1.4 (as in the Blandford 
& Konigl (1979) and MFF models), the radio luminosity is expected to dominate the X-rays below L x ~ 4 X 10- 5 L Edd (about 0.005 
Crab). Possible detection of the nearby SXT A 0620-00 at the predicted radio level would be capable to probe such a statement at very 
low X-ray level, requiring a wholly new accretion regime on to stellar BHs. 



5 BEYOND THE HARD-TO-SOFT STATE 
TRANSITION 

5.1 Discrete ejections 

So far, we only have focused on low/hard state BHs, that 
is, binary systems characterised by a quasi-steady state of 
stable accretion and whose X-ray spectrum is dominated by 
a hard power law. 

In the following we will add to our sample radio and X-ray 
fluxes from transient BHCs during their episodic outbursts 
associated with discrete ejection events corresponding to op- 
tically thin radio emission. 

All the available data, comprising of soft X-ray and radio 



peak fluxes with references, have been reported by Fender & 
Kuulkcrs (2001). Two points we include in this Section come 
from systems which have also been shown in Figure 2 while 
in the hard state, namely V 404 Cygni and GRO J0422+32. 
Sources under consideration, as well as their fundamental 
physical parameters, are listed in Table 3. Simultaneous 
radio:X-ray peak fluxes, scaled to a distance of 1 kpc are 
plotted in Figure 6 with filled squares. 

Radio data are based on peak observed flux density at fre- 
quency of 5 GHz. Where measurements at 5 GHz were not 
available, a spectral index of a = —0.5 was assumed in order 
to estimate the 5 GHz flux based on observations performed 
at different frequencies (see Fender & Kuulkers 2001 for de- 
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Table 3. Transient (plus 4 peculiar) sources whose simultane- 
ous radio:X-ray peak fluxes associated with discrete ejections are 
plotted in Figure 6 (data from Fender & Kuulkers 2001). 



Source 


Distance 


(10 22 CTO~ 2 ) 


1A 0620-00 


1 (1) 


0.2 (a) 


GS 2000+25 


2 (2) 


0.2 (b) 


GS 1124-68 


4.9 (3) 


0.5 (c) 


GRO J1655-40 


3(4) 


0.7 (d) 


GRS 1716-249 


2.4 (5) 


0.4 (e) 


GRS 1739-278 


8.5 (6) 


2.6 (f) 


4U 1630-472 


10 (7) 


9.4 (g) 


XTE J2012+381 


10 (8) 


1.3 (h) 


XTE J1748-288 


8 (9) 


7.5 (i) 


XTE J1859+226 


7.6 (10) 


0.3 (1) 


XTE J0421+560 


5 (11) 


0.1 (m) 


GRS 1915+105 


11 (12) 


7.0 (n) 


Cygnus X-3 


9 (13) 


2.7 (o) 


LMC X-l 
LMC X-3 


55 (14) 
55 (14) 


0.8 (14) 
0.1 (14) 



References : 1: Shahbaz et al. 1994; 2: Callanan et al. 1996; 
3: Shahbaz et al. 1997; 4: Kubota et al. 2001; 5: Delia Valle et 
al. 1994; 6: Marti et al. 1997 ; 7: Parmar et al. 1986; 8: 
Campana et al. 2002; 9: Kotani et al. 2000; 10: Hynes et 
al. 2002; 11: Robinson et al. 2002; 12: Fender et al. 1999a; 13: 
Predehl et al. 2000; 14: Haardt et al. 2001. 

a: Kong et al. 2002; b: Rutledge et al. 1999; c: Ebisawa et 
al. 1994; d: Ueda et al. 1998; e: Tanaka 1993 f: Greiner et 
al. 1996; g: Tomsick & Kaaret 2000; h: Hynes et al. 1999; i: 
Miller et al. 2001; 1: Markwardt et al. 1999; m: Parmar et 
al. 2000; n: Klcin-Wolt et al. 2002; o: Tcrasawa & Nakamura 
1995. 

tails). 

Most of the X-ray data are ASM detections, with the 
only exceptions of GRO J1655-40 and GRO J0422+32, 
whose outbursts have been detected by either BATSE or 
GRANAT. For clarity, no error bars are plotted. 

5.2 Other BHCs: persistent soft state and 
'extreme' sources 

In between the jet quenching and the discrete ejections from 
transient sources, there are of course binary systems display- 
ing a persistent soft spectrum, i.e. whose emission is domi- 
nated by disc blackbody photons. LMC X-l and LMC X-3, 
in the Large Magellanic Clouds, are the only BHCs always 
observed while the soft state (actually Wilms et al. 2001; 
Boyd & Smale 2000 and Homan et al. 2000 reported signs 
that LMC X-3 entered in a hard state). For both these 
sources, the presence of a black hole is quite well established, 
with a most likely mass of 9Mq and 6Mq, respectively for 
LMC X-3 and LMC X-l. By means of radio and X-ray ob- 
servations performed in 1997, we can place the two sources 
in the S ra dio vs. Sx plane, in order to verify the amount of 
radio power from soft state BHs, which is expected to be 
well below the hard state correlation. LMC X-3 X-ray flux 
densities have been derived from Haardt et al. (2001), while 
LMC X-l from Gierlinski et al. 2001. Radio upper limits for 



both sources are taken from Fender, Southwell & Tzioumis 
(1998). The corresponding values, scaled to a distance of 
1 kpc and corrected for absorption according to equation 
2, are plotted in Figure 6 with filled circles (at about 150 
Crab:<4540 mjy and 45 Crab:<360 mjy, LMC X-l and 
LMC X-3 respectively). As expected, both points lie below 
the hard state relation extended to such high X-ray ener- 
gies. Although we are only reporting upper limits, we can 
assert that fluxes from LMC X-l and LMC X-3 do not dis- 
agree with our previous finding. 

For completeness, simultaneous radio:X-ray detections of 
two extreme sources, namely GRS 1915+105 (open poly- 
gons) and Cygnus X-3 (open circles), have been included, 
corresponding to the two big 'clouds' at the top right re- 
gion of Figure 6. Both systems are traditionally considered 
'exotic' due to their timing and spectral behaviours which 
do not fully resemble any 'standard' picture generally ac- 
cetpted for BHCs; for instance, both these sources display 
either optically thin or thick radio spectra. For extensive re- 
views, see Bonnet-Bidaud & Chardin (1988; Cygnus X-3) 
and Belloni et al. (2000; GRS 1915+105). Here we note that, 
despite its unusually high luminosity, detections of GRS 
1915+105 in the so called plateau state (Belloni et al. 2000) - 
which appears to share similar properties with the canonical 
low/hard state - still seem to belong to the 0.7 correlation 
extended up to super-Eddington regime. 

Cygnus X-3 is the strongest observed persistent radio- 
emitting BH binary and is embedded in a dense stellar 
wind from the companion Wolf-Rayet star (van Kerkwijk 
et al. 1992; Fender, Hanson & Pooley 1999), which makes 
difficult to isolate the compact object high energy spectrum. 
The high energy emission from the vicinity of the compact 
object in Cygnus X-3 is likely to be hidden by a dense stellar 
environment surrounding the source; as a consequence, the 
intrinsic X-ray luminosity might be higher than inferred, 
pushing the dataset closer to the 0.7 relation. Moreover, 
there still remains uncertainty about the nature of the ac- 
cretor in this system. The neutron star hypothesis can not 
be ruled out with confidence. 

It is interesting to note that, while the jets from GRS 
1915+105 are at 60-70° (Fender et al. 1999a), those of 
Cygnus X-3 appear to be close to the line of sight ( < 14°, 
Mioduszewski et al. 2001), supporting the previous hypoth- 
esis, in which higher-than-average normalisation factors 
would correspond to higher Doppler factors. In addition, the 
behaviour of Cygnus X-3, with the apparent turnover in the 
radio power, is very similar to that of Cygnus X-l, except 
for the flaring behaviour following the jet quenching (note 
that points corresponding to this flaring show characteristic 
optically thin radio spectra; on the contrary, pre-flaring de- 
tections display 'flat' radio spectra, i.e. a different physical 
origin). 



6 PREDICTING RADIO FLUXES AT LOW 
QUIESCENT LUMINOSITIES 

As the same correlation appears to be maintained over many 
years and for different sources (like for instance Cygnus X-l, 
GX 339-4, V 404 Cyg), we can estimate the level of radio 
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Table 4. Five X-ray transients have been monitored in X-rays by different telescopes in different energy ranges (indicated below the 
table) during quiescence. Here we report the values for the maximum and the minimum luminosity, with relative references and energy 
ranges. For these values the corresponding predicted radio flux densities have been calculated based upon the radio:X-ray correlation we 
have found, that is S ra di = [223 X (Sx,ikpc/Crab) 0.7]/ (D /kpc) 2 m.]y. 



Source 


X-ray luminosity 


Distance 


Flux density 


Predicted radio 




(10 32 erg/sec) 


(kpc) 


at 1 kpc (10~ 6 Crab) 


flux density (/iJy) 


1A 0620-00 


0.02 1 - 0.04 2 (a,b) 


1 


1-5 


18-43 


GRO J1655-40 


0.2 3 - 3 4 (a,c) 


3 


6-82 


6-34 


XTE J1550-564 


< 5 5 (d) 


4 


< 173 


< 32 


GRO J0422+32 


0.08 4 (c) 


2.4 


~ 2 


~4 


GS 2000+25 


0.02 4 (c) 


2 


~ 0.5 


~2 



References : 1: 0.4-2.4 keV; 2: 0.4-1.4 keV; 3: 0.3-7 keV; 4: 0.5-10 keV; 5: 0.5-7 keV. 

a: Kong et al. 2002; b: Narayan et al. 1996; c: Asai et al. 1998; d: Tomsick et al. 2001; e: Garcia et al. 2001. 



emission from a hard state BH by measuring its X-ray flux 
alone. This is particularly interesting in the case of black 
hole X-ray transients, whose inferred accretion rate during 
quiescence may be very small. 

Kong et al. (2002) present Chandra observations of three 
BH transients during quiescence for which no simultaneous 
radio detection is available to date, namely A 0620-00, GRO 
J1655-40 and XTE J1550-564. In order to check for possi- 
ble spectral variability, they also report results from previous 
X-ray observations carried out by different telescopes, such 
as ROSAT, ASCA and BeppoSAX. According to Tomsick 
et al. 2001, the lowest luminosity measured for XTE J1550- 
564 with Chandra (5 x 10 32 erg/sec, for a distance of 4 kpc) 
should however be considered only as an upper limit on the 
quiescent luminosity of the system. Chandra detections of 
other two transient sources, namely GRO J0422+32 and GS 
2000+25, are reported by Garcia et al. 2001. In Table 4 we 
list for each of the five sources the maximum and the mini- 
mum measured X-ray luminosity in quiescence, the inferred 
distance and the predicted radio flux density (in pJy) based 
on our best-fit equation, that is S ra dio = 223 x (Sx) +0 ' 7 ■ 
Given a spread of 156 over 223 in the normalisation factor, 
the predicted values must be considered reliable within one 
order of magnitude. 

6.1 A 0620-00: the ideal candidate 

The Soft X-ray Transient (SXT) A 0620-00 was discov- 
ered in outburst in 1975 August (Elvis et al. 1975), while 
the associated radio source was at a level of 200-300 mjy 
(Owen et al. 1976) during the onset of the outburst. Six 
years later the source was detected with the VLA at level of 
249±79^Jy (Geldzahler 1983); additional VLA observations 
in 1986 (see McClintock et al. 1995) yielded an upper limit 
of 140/xJy, clearly indicating a decline in radio power (see 
Figure 6 where A 0620-00 detection/upper limit/predicted 
radio fluxes are marked with triangles). The 1981 detection 
might actually be associated to radio lobes resulting from 
the interaction of a relativistic-decelerating jet with the in- 
terstellar medium, as observed in the case of XTE J1550-564 
about four years after the ejection of plasma from near the 
BH (Corbel et al. 2002). 

Due to its relative proximity, A 0620-00 is the most suitable 
candidate to probe if our empirical radio:X-ray relation does 



hold down to low quiescent X-ray luminosity (~ 2 x 10 30 
erg/sec at 1 kpc, i.e. about W~ 8 L E dd', Garcia et al. 2001). 
In other words, if A 0620-00 was detected at the predicted 
radio level (a few tens of /iJy, see Figure 6, bottom left cor- 
ner), it would confirm that the mechanisms at the origin of 
radio and X-ray emission are correlated, if not even partly 
coincident, over more than six orders of magnitude in X-ray 
luminosity. 

Moreover, if the radio:X-ray correlation were confirmed at 
very low X-ray luminosities (below 10 4 L E dd ) it would 
strongly constrain the overall theory of accretion in quies- 
cence. We direct the reader to the discussion for further 
comments. 



7 DISCUSSION 

The presence of a coupling between radio and X-ray emis- 
sion in the low/hard state of black hole binaries obviously re- 
quires a theoretical interpretation that relates somehow the 
powering/quenching mechanism(s) of the jet to the overall 
accretion pattern. Zdziarski et al. (2003) ascribe the correla- 
tion to a correspondence between the level of X-ray emission 
and the rate of ejection of radio-emitting blobs forming a 
compact jet. In this picture, there still remains the ques- 
tion of the condition for jet suppression. Meier (2001) inter- 
prets the steady-jet/hard-X-ray state association as strong 
evidence for magnetohydrodynamic (MHD) jet formation, 
where the most powerful jets are the product of accretion 
flows characterised by large scale height. The simulations 
show in fact that the jet is confined by the toroidal compo- 
nent of the magnetic field lines, coiled due to the disc differ- 
ential rotation: the bigger the disc scale-height, the stronger 
the field. The power of the jet naturally decreases (at least 
100 times weaker) in the soft/high state, associated with a 
standard geometrically thin disc (Shakura & Sunyaev 1973). 
To date, two broad classes of geometries have been pro- 
posed for explaining X-ray emission from low/hard state 
BHs. The more classical picture is that of a hot, homo- 
geneous, optically thin corona of (quasi-)thermal electrons 
which inverse Compton scatter 'seed' photons coming from 
the underlying accretion disc (Shapiro, Lightman & Eard- 
ley 1973 and similar later solutions) . The alternative is that 
of a jet -synchrotron model, in which, under reasonable as- 



12 E. Gallo, R. P. Fender & G. G. Pooley 



sumption on the input power of the jet and the location 
of the first acceleration zone, optically thin synchtrotron 
emission can dominate the X-ray spectrum, traditionally 
thought to be a product of inverse Compton process only 
(see Markoff, Falcke & Fender 2001 for a detailed description 
of the model; hereafter MFF). This model however predicts 
self-synchrotron Compton upscattering in the jet for some 
scenarios. Remarkably, the MFF model is able to reproduce 
the observed slope of the radio:X-ray correlation analyti- 
cally (Markoff et al. 2003a), as a function of the measured 
X-ray and radio-infrared spectral indices. 

A revised version of the classical Comptonisation model for 
the hard state has been proposed by Beloborodov (1999). 
In this case the hot coronal plasma is powered by magnetic 
field line reconnections (Galeev, Rosner & Vaiana 1979) and 
confined within several active flares with mildly relativistic 
bulk velocities, inferred by the relative weakness of the re- 
flection component. Due to aberration effects in fact, the 
amount of X-rays as seen by the reflecting disc turns out to 
be reduced by a factor consistent with f3x — 0.3. 
Maccarone (2003) draws a similar conclusion on a different 
ground: he has tabulated all the the available measuremets 
of X-ray luminosities at the time of the soft -to-hard state 
transition, for both BHC and neutron star systems. The 
resultant variance in state transition luminosity is also con- 
sistent with coming from material with f3x — 0.2; therefore 
Fx — 1 in both cases. 

Following the approach of Section 4.2.2, this immediately 
implies a stringent upper limit on the beaming of radio 
emission as well, that is r radio £ 2 (see Figure 5). This is 
almost certainly significantly less relativistic than the jets 
produced during transients outbursts of sources such as GRS 
1915+105 (Mirabel & Rodriguez 1994; Fender et al. 1999a) 
and GRO J1655-40 (Hjellming & Rupen 1995; Harmon et 
al. 1995; see also Fender 2003). 

Therefore, if mildly relativistic beaming characterises the 
low/hard state, a mechanism(s) must exist which both 
switches the jet off - high/soft state - and produces a faster 
jet - discrete ejections - above ~ W~ 2 Lsdd , where the 
hard-to-soft state transition occurs. 

It is interesting to mention that a few percent of the Ed- 
dington rate is also close to the regime at which Ghisellini 
& Celotti (2001) have identified the transition line between 
FRI and FR II radio-galaxies, the former class being associ- 
ated with slower kpc-scale jets than the latter (see Begelman 
1982; Bicknell 1984; Laing 1993). This leads to a kind of cor- 
respondence between 'extreme sources' like GRS 1915+105, 
or Cyg X-3, and FRII, characterised by quite high bulk 
Lorcntz factors, while 'canonical' hard state sources would 
be associated with FR I, with relatively low outflow bulk 
velocities. 



7.1 Jet— dominated 'quiescence' 

A main task of the models remains that of reproducing the 
observational behaviour of accreting stellar BHs at a vari- 
ety of accretion rates. An interesting prediction of this work 
concerns the relative power of the jet, with respect to the 
overall accretion power, at low quiescent luminosities. 
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Figure 7. Broadband spectra of Cygnus X-l, V 404 Cygni and 
GX 339-4: for low/hard state (LHS) black holes, not only the 
behaviours in the S ra( n / Sx plane look very similar, but even 
the shape of their energy spectrum, from radio wavelengths up 
to X-rays, suggesting cither a common origin or coupling of the 
basic emission processes. 



In this regime, the BH spectral energy distribution appears 
very similar to that of canonical hard state sources (see Fig- 
ure 7), although the X-ray emission in quiescence is not well 
reproduced by the standard accretion-corona model, requir- 
ing a much lower radiative efficiency. 

Narayan et al. (1996, 1997, 2001) showed that an Advec- 
tion Dominated Accretion Flow (ADAF), in which the en- 
ergy released by viscous torques is 'stored' into the flow 
rather than radiated away, can adequately model the avail- 
able observations at high energies. However, on the other 
extreme of the spectrum, in the radio band, the existence 
of a jet seriously weakens such a solution, requiring a sig- 
nificantly different physics to model the observed spectrum. 
Jet-powered radio emission is in fact several orders of mag- 
nitude brighter than expected extrapolating ADAF spectra 
down to radio band; moreover, the standard ADAF picture 
predicts highly-inverted radio spectra, instead of the ob- 
served flat ones. 

An alternative scenario for low-luminosity stellar black holes 
has been proposed by Merloni & Fabian (2002). They show 
that a coronal-outflow dominated accretion disc, in which 
the fraction of the accretion power released in the corona 
increases as the accretion rate decreases, would be an ideal 
site for jet-launching, both MHD and thermally driven. 

Based on the radio:X-ray correlation, the ratio between the 
observed radio and X-ray luminosities scales as: 



-"radio / T \ —0.3 

~Lx~ X 



(8) 
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This already implies that, as the X-ray luminosity decreases, 
the radiative jet power will become more and more impor- 
tant with respect to the X-rays. Moreover, because of the 
self-absorption effects, it has been shown that there will 
be not a linear relationship between the radio luminosity, 
Lradio, and the total jet power, Lj et , for any optically thick 
jet model which can explain the fiat radio spectrum observed 
during the hard state. Blandford & Konigl (1979), Falcke & 
Biermann (1996), and Markoff, Falcke & Fender (2001) ob- 
tain in fact L ra dio oc (Lj et ) 1A . 

If so, and in general for any relationship of the form: 

Lradio OC (Lj et ) X (9) 

equation 8 implies the following scaling for the fractional jet 
power: 

^ oc [L x f^ (10) 

Hence, for any x > 0.7, there exists an X-ray luminosity 
below which the jet will be the dominant output channel for 
the accretion power. 

If x = 1.4, by re-scaling the numbers to XTE J1118+480 - 
emitting at ~ 10~ L E dd, and whose fractional jet power is at 
least 20% (Fender et al. 2001) - one obtains that Lj et Lx 
for L x ~ 4 x 10~~'L E dd (i-e. below 0.005 Crab - scaled - 
see Figure 6). Note that observations of GX 339-4 and V 
404 Cygni in quiescence already cover this regime, meaning 
that, if x — 1.4, both these sources actually are in a jet- 
dominated state. 

This would be a wholly new accretion regime for X-ray 
binaries, requiring significant modification of existing (e.g. 
ADAF) models. In addition it would indicate that the 
overwhelming majority of 'known' stellar-mass black holes, 
which are currently in quiescence, are in fact feeding back 
most of their accretion energy into the interstellar medium 
in the form of the kinetic energy of the jets and are accret- 
ing at rather higher levels than derived based only on their 
X-ray luminosity. 

If we do establish that accretion is taking place in quies- 
cence, for instance throught the detection of A 0620-00, and 
is furthermore channelling most of its power into jet forma- 
tion, then the arguments for observational evidence for black 
hole event horizons based upon a comparison of quiescent 
X-ray luminosities of black hole and neutron star binaries 
(e.g. Garcia et al. 2001) will need to be re-examined. In 
fact, assuming that Lj et cc (Lx)" 0,7/ ^' holds for both neu- 
tron stars and black holes, then the observed difference in 
'radio loudness' between black hole and neutron star bina- 
ries (Fender & Kuulkers 2001) might be enough on its own 
to explain the discrepancy, and it may be that the event 
horizon plays no part. This is explored further in Fender, 
Gallo & Jonker (in preparation). 



8 SUMMARY 

In this paper we provide observational evidence for a broad 
empirical relation between radio and X-ray emission in 
Galactic black hole binaries during their spectrally hard 
state. The main points established throught this work can 
be summarised as follows: 

• In low/hard state BHCs the observed radio and X-ray 



fluxes are correlated over more than three orders of magni- 
tude in accretion rate, with a spread in radio power of about 
one order of magnitude. 

• Even at accretion rates as low as 10~ 5 Eddington a 
powerful jet is being formed; no lower limit to the relation 
has been found. 

• V 404 Cygni is the second source to display S ra dio oc 
(Sx) 0,7 , from quiescence up close to the hard-to-soft state 
transition. 

• Assuming 0.7 as a universal slope for the low/hard 
state, and under the hypotheses of a) common disc-jet cou- 
pling and b) isotropic X-ray emission, the measured spread 
in radio flux can be interpreted in terms of a distribution 
in Doppler factors. Monte Carlo simulations show that the 
observed scatter is consistent with relatively low beaming 
(Tradio ~ 2) outflows in the low/hard state, unlike those in 
transient outbursts. 

• When the combination of radio and X-ray beaming is 
taken into account, the range of possible bulk velocities in 
the jet significantly broadens, allowing the X-ray emitting 
material to be relativistic for almost any value of Y ra dio, 
but implying strong X-ray selection effects. In this case an 
independent estimation on Fx is needed to limit T 'radio- Un- 
related works (Beloborodov 1999; Maccarone 2003) impose 
stringent constraints on the bulk velocity of the X-ray emit- 
ting material, leading to the conclusion of relatively low ra- 
dio beaming {Y ra dio ~ 2) in the hard state. 

• Close to the hard-to-soft state transition the jet 
switches off, probably in all sources. The X-ray luminosity 
at which the radio quenching occurs might positively corre- 
late with the BH mass, being consistent with taking place at 
a constant fraction of the Eddington rate. It is worth men- 
tioning that a similar fraction of L E dd has been identified as 
a dividing line between FR I and FR II radio-galaxies, that 
is between supermassive BHs producing mildly and highly 
relativistic jets respectively (Ghisellini &: Celotti 2001). 

• Since the correlation appears to be maintained over 
many years and for different sources, this leads to the pos- 
sibility of predicting the level of radio emission from a hard 
state and/or quiescent BH by measuring its X-ray flux. 

• If the radio luminosity scales as the total jet power 
raised to x, with x > 0.7, this implies the existence of an X- 
ray luminosity below which the most of the accretion power 
will be channelled into the jet rather than in the X-rays. If 
x = 1.4 (e.g. Blandford & Konigl 1979, Falcke & Biermann 
1996; MFF), then below L x ~4x 10~ 5 L Edd the jet is 
expected to dominate. 

This work provides evidence for a physical coupling between 
radio and hard X-ray emitting ouflows from accreting stellar 
BHs. A key, still unresolved issue concerns the modelling of 
the transition between X-ray states in a self consistent way, 
which could possibly account for both the jet suppression, 
when the disc dominates, and the transition from mildly to 
highly relativistic jets, as in case of transient outbursts. In- 
cluding the formation of jets in the overall energetics and 
dynamics of the accretion process at a variety of X-ray lu- 
minosities has undoubtedly become of primary importance 
to address, especially based on mounting evidence for the 
jet power to be a significant fraction, if not the dominant 
output channel, of the total accretion power. 
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